TWO-STREAM INSTABILITY OF COUNTER-ROTATING 

GALAXIES 



R.V.E. Lovelace, K.P. Jore, and M.P. Haynes 
Department of Astronomy, Cornell University, Ithaca, NY 14-853-6801 



ABSTRACT 

The present study of the two-stream instability in stellar disks with counter-rotating 
components of stars and/or gas is stimulated by recently discovered counter-rotating spi- 
ral and SO galaxies. Strong linear two-stream instability of tightly-wrapped spiral waves 
is found for one and two-armed waves with the pattern angular speed of the unstable 
waves always intermediate between the angular speed of the co-rotating matter (+0) 
and that of the counter-rotating matter (— fl). The instability arises from the interaction 
of positive and negative energy modes in the co- and counter-rotating components. The 
unstable waves are in general convective - they move in radius and radial wavenumber 
space - with the result that amplification of the advected wave is more important than 
the local growth rate. For a galaxy of co-rotating stars and counter-rotating stars of 
mass-fraction < |, or of counter- rotating gas of mass- fraction £ g < |, the largest 
amplification is usually for the one-armed leading waves (with respect to the co-rotating 
stars). For the case of both counter-rotating stars and gas, the largest amplifications 
are for £* + £ g ~ h, also for one-armed leading waves. The two-armed trailing waves 
usually have smaller amplifications. The growth rates and amplifications all decrease 
as the velocity spreads of the stars and/or gas increase. It is suggested that the spiral 
waves can provide an effective viscosity for the gas causing its accretion. 
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1. INTRODUCTION 

A surprise of recent high spectral resolution stud- 
ies of normal galaxies is the occurrence of counter- 
rotating gas and/or stars in galaxies of all morpho- 
logical types, ellipticals, spirals, and irregulars (see re- 
view by Gallctta 1996). In the ellipticals, the counter- 
rotating component is usually in the nuclear core and 
may result from merging of galaxies with opposite an- 
gular momentum. In contrast, in a number of spiral 
and ,50 galaxies, counter-rotating disks of gas and/or 
stars have been found to co-exist with the primary 
disk out to large distances (10 — 20 kpc). Examples 
include NGC 4550 (Rubin, Graham, & Kenney 1992; 
Rix et al. 1992), NGC 4826 (Braun ct al. 1992), 
NGC 7217 (Merrifield & Kuijken 1994), NGC 4546 
(Sage & Galetta 1994), NGC 3626 (Ciri et al. 1995), 
NGC 3593 (Bertola et al. 1996), and NGC 4138 (Jore, 
Brocils, & Haynes 1996). Table 1 summarizes data on 
these galaxies. These galaxies are early type and have 
ring(s) or dust lanes or other morphological peculiar- 
ities. Related cases include galaxies with infalling HI 
gas, for example, NGC 4254 which has an m — 1 spi- 
ral arm (Phookun, Vogel, & Mundy 1993), and NGC 
628 which has a warped and distorted outer disk with 
high- velocity clouds (Kamphius & Briggs 1992). 

It is not likely that the large-scale counter-rotating 
components result from mergers of similar mass galax- 
ies of opposite angular momenta because of the large 
vertical thickening observed in simulation studies of 
such mergers (Barnes 1992). Formation of disk galax- 
ies by accretion has been studied by Ryden and Gunn 
(1987) and Ryden (1988). Quinn and Binney (1992) 
have shown that the more recently accreted mat- 
ter at large distances has a spin anticorrelated with 
the spin of the central mattter. Thakar and Ryden 
(1996) discuss different possibilities for the forma- 
tion of counter-rotating galaxies, (1) that the counter- 
rotating matter may come from the merger of an op- 
positely rotating gas rich dwarf galaxy with an exist- 
ing spiral, and (2) that the accretion of gas may occur 
over the lifetime of a galaxy with the more recently 
accreted gas counter-rotating. 

An important open theoretical question is: What is 
the interaction between the co- and counter-rotating 
components observed on large-scales in spiral and 5*0 
galaxies, and docs this interaction facilitate accre- 
tion of the counter-rotating component? The com- 
ponent without neutral or ionized gas rotating in the 
same direction is referred to as the co-rotating com- 



ponent. Clearly, a strong interaction will occur be- 
tween co- and counter-rotating gas with the result 
that these components are likely to be spatially sep- 
arated (see however Lovelace and Chou 1996). In 
this work we assume that any co-rotating gas has 
been swept-out from the spatial regions (radii) occu- 
pied by counter-rotating gas. The interaction between 
co- and counter-rotating stars by dynamical friction 
is negligible because of the small densities and large 
relative velocities. Also, the direct interaction of co- 
rotating stars with counter-rotating gas is negligible 
because of the small cross-section of a star. Accre- 
tion occurs in the counter-rotating galaxy model of 
Thakar and Ryden (1996) due to viscosity of the gas 
which is treated using sticky-particle TV-body simu- 
lations. However, the viscosity used is not derived 
from physical principles. Here, we analyze the influ- 
ence of the two-stream instability in generating spiral 
waves in counter-steaming flat galaxies and suggest 
that these waves may give rise to an effective viscos- 
ity for counter-rotating gas. 

Theoretical interest in galaxies with both co- and 
counter-rotating stars was stimulated by the early sta- 
bility analysis of Kalnajs (1977) who found that the 
counter-rotating stars have a stabilizing influence on 
the bar forming (m = 2) mode in models without 
a massive halo. Studies of counter-streaming stabil- 
ity of idealized self-gravitating systems were made 
by a number of authors, by Bisnovatyi-Kogan ct 
al. (1969) and Bisnovatyi-Kogan (1973) for counter- 
rotating stars in a finite radius cylinder, by Araki 
(1987) for equal populations of counter-rotating stars 
in finite-radius rigidly rotating Kalnajs disks, and by 
Lyndcn-Bcll (1967) and Araki (1987) for a homo- 
geneous system of counter-streaming stars following 
Jeans' neglect of the background potential. Pertur- 
bations of the form f(r)exp(im<p — iuit) (in cylindri- 
cal coordinates (r, <f), z) with m = 0, 1, ..) were found 
to be linearly unstable for some conditions for the 
cylinder for m = 2, whereas for the Kalnajs disks 
the m = 1 mode appears as the main two-stream 
instability. For homogeneous counter-streaming sys- 
tems with Maxwellian distribution functions, Lynden- 
Bcll and Araki found that the Jeans instability al- 
ways dominates the two-stream instability. N-body 
computer simulations of stellar disks by Zang and 
Hohl (1978), Sellwood and Merrit (1994), and Howard 
ct al. (1996) support the idea of Kalnajs (1977) 
that counter-rotating stars have a stabilizing influ- 
ence on the bar forming mode. However, the counter- 
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streaming is found to enhance the growth of the one- 
armed (m = 1) mode. 

Section 2 of this work discuss the basic equations 
for the small-amplitude, tightly-wrapped spiral waves 
in a flat galaxy. Section 3 treats the stability of 
the spiral waves for different cases of a disk with co- 
rotating stars and counter- rotating gas and/or stars. 
Section 4 shows that the two-stream instability found 
in Sec. 3 arises from the interaction of positive and 
negative energy modes associated with the co- and 
counter-rotating matter. Section 4 discusses possible 
non-linear effects of the instabilities including gas ac- 
cretion. Section 5 gives conclusions of this work. 

2. THEORY OF SPIRAL WAVES 

We first give a brief summary of the linear WKB 
theory of tightly-wrapped spiral waves in a single 
component galaxy of stars rotating with angular rate 
f2(r) > (sec for example Binney & Tremaine 1987; 
hereafter denoted BT; or Palmer 1994). We use an 
inertial, cylindrical (r, <j), z) coordinate system and as- 
sume a thin disk galaxy as indicated in Figure 1. In 
the midplane of the galaxy the perturbation of the 
gravitational potential is 

<5$(r, <f>, z = 0,t) = C cxp [ i J dr' k r (r')+im<p—iu)t\ , 

(1) 

where the radial wavenumber k r satisfies \ k r r\ ^> 1 for 
a tightly wrapped wave, C(r) a slowly varying func- 
tion of r, m = ±1, ±2, .. is the number of spiral arms, 
and oj is the angular frequency of the wave. Only 
m > need be considered because 5$* is a valid so- 
lution if (5$ is. Then, k r > (< 0) corresponds to a 
trailing (leading) spiral wave. The thin disk assump- 
tion requires \k r h\ <C 1, where h is disk half-thickness. 

A well-known calculation (see BT or Palmer 1994) 
of the linear perturbation of the dynamical equations 
leads to the dispersion relation 

= e(u),k r ) = l+V4cj,k r ) , (2) 

where 

_ 2|fc r |exp(-X.) v h(X«) 

s=iuj-mQ)/K, k 2 = \-^-{r 4 n 2 ) , 

n 2 

k r — k r J k cr it , kcrit — ~Z 7^ ' 



X, = = 0.28568(QX) 2 , 

= K<T r 

^* ~ 3.3583GS tot ' 

Here, e(u>, k r ) has the role of a dielectric function 
for the disk, and V* the polarization function (the 
♦-subscript denotes stars); s is a dimcnsionlcss fre- 
quency; n(r) is the radial epicyclic frequency of a star; 
k cr u is a characteristic wavenumber, and k r is the di- 
mensionless radial wavenumber; T,tot(r) is the total 
surface mass density of the disk; Q*(r) is Toomre's 
stability parameter; oy(r) is the radial dispersion of 
the star velocities (for a stellar distribution function 
/* cx cxp[— v 2 /(2a 2 )]); and is the usual modified 
Bcsscl function of order The condition for axisym- 
metric (m = 0) stability is (Q* > 1 (Toomre 1964). 
The influence of halo matter enters through fl(r). 

Figure 2 shows wavenumber (k r ) - frequency (s) 
plots for two values of Q» obtained from equation 
(2). The corresponding plot in BT (Figure 6-14a) 
and in Palmer (1994, Figure 12.2) is incomplete be- 
cause the higher order branches, labeled / = 2,3,.., 
have been omitted. Further, the statement in BT 
that "a stellar disk has no pressure forces and there- 
fore cannot support waves with \s\ > 1" (page 367) 
is incorrect. The branches I — 2,3,.. are the ana- 
logues of the well-known Bernstein modes which prop- 
agate across a uniform magnetic field in a collisionless 
plasma (sec for example Krall & Trivelpiece 1973). 
Here, the coriolis force is the analogue of the Lorentz 
force in a magnetized plasma. Note in particular 
that the Linblad resonances correspond to s = ±1 
or flp = Q ± Ik/tu, where f2 p = cu/m is the pattern 
angular speed, m = 1, 2, .., and I = 1,2, ... 

For later use, we also give the dispersion relation 
for a single component gaseous disk with rotation rate 
O(r) > 0, 

Q = e{uo,k r ) = l+V g {uo,k r ) , (3) 

where 

\K\ 

9 S 2 - 1 - {k r C s /K) 2 ' 

where c s is the sound speed in the gas, and where 
the other quantities are the same as in equation (2). 
Equation (3) gives the explicit dependence s 2 = 1 — 
\k r \ + (k r c s /n) 2 . Here, the Toomre (1964) condition 
for axisymmetric stability is c s > n/(2k cr it) or Q g = 
Kc s /(7rGS tot ) > 1. 
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Note that f2 2 , k 2 , and k = |/t| are the same for 
both co- and counter-rotating disk components. Note 
also that S tot is the total - gas plus stellar - mass 
density in all components. In that P* and V g are 
even functions of k r , we have u>(k r ) = uj{—k r ). 

3. TWO-STREAM INSTABILITY 

In the following subsections we consider the two- 
stream instability in a number of different limits. 
Specifically, we consider co-rotating stars (+fi) and 
in subsection 

3.1 counter-rotating gas (— fi), 

3.2 a small mass-fraction of counter-rotating gas, 

3.3 a large mass-fraction of counter-rotating gas, 

3.4 a small mass-fraction of counter-rotating stars, 

3.5 a large mass- fraction of counter-rotating stars, 
and 

3.6 arbitrary fractions of counter-rotating gas and 
stars. 

3.1. Co-Rotating Stars/Counter- Rotating Gas 

The dispersion relation for a two component galaxy 
consisting of co-rotating stars (+fi) and counter- 
rotating gas (— fi) can be written down immediately 
using equations (2) and (3), 

= e(w,k r ) = l + (l-^ g )V*(s,k r )+^ g V g (s+w,k r ) , 

(4) 

where w = 2mQ/n , and where 

is the fraction of the disk surface mass density in gas. 
To understand equation (4) it is useful to introduce 
the notion of stellar modes which obey 

= l + (l-tg)V,(8,kr) , (5a) 

and gas modes which obey 

= l + £ g V g {s + w,k r ) . (56) 

Equation (5a) gives a family of curves s — s*(fc r ) for 
the star modes, while (5b) gives curves s — s g (k r ) for 
the gas modes. This is shown in Figure 3 for m = 2. 

The approximation of equations (5) breaks down 
near the isolated points in the (k r , s) plane where the 
star and gas modes cross, s*(fc r ) = s g (k r ). At these 
points there is a strong resonant interaction between 
the modes. 



3.2. Co-Rotating Stars/Low-Mass Counter- 
Rotating Gas 

For ( 9 < 1, an approximate solution of the dis- 
persion relation (4) near a mode-crossing point can 
easily be obtained as follows. To zeroth order in 
£ g , equation (4) is satisfied by (k r ,s a ) obeying = 
1 + V*(s , k r ). For small £ s , the gas response is large 
if e = [(.s + w) 2 - 1 - (k r c s /n) 2 ]/[2{s + w)} is small. 
To first order in £ s , the solution of equation (4) can 
thus be written as s = s + 5s with \5s\ <C \s a \. Taylor 
expanding about s gives 

5s = — = -e±(e 2 + 5s 2 )\ (6) 

where 

t I h I 

S S 2 = ^ IglJ 

°- 2(8 + w)(&P./d8 ) ' 

Instability, u>i = Im(uj) > 0, occurs if Ss 2 < and 
e 2 < \Ss 2 \. (The damping roots with Wi < are 
ignored here and subsequently.) From equation (2), 
dV*jds = — s a \..\. Thus for instability s (s + w) < 
0. Equivalcntly, there is instability if the frequency of 
the mode crossing, w = ns Q + mCl, obeys 

-n < — < n , (7) 

m 

where cj /m is the pattern velocity of the perturba- 
tion. Later, in Sec. 4, we give a more general treat- 
ment of instability in counter rotating galaxies which 
also leads to equation (7). From equation (6), note 

that the maximum growth rate, max(iVi) = k\5s \, 

i 

occurs for e — and scales as Q . 

Figure 4 shows the behavior of the mode cross- 
ings for an unstable and a stable crossing of Fig- 
ure 3 obtained from equation (4). For £ g -c 1 and 
Csjor = 1, there are no unstable crossings for m = 1, 
while for m = 3 there are unstable crossings near 
k r rts 1.04, 2.9, and 4.24. In all cases the unstable 
crossings obey equation (7). 

Consider the consequence of the wave growth. For 
this it is useful to examine the evolution of a wave 
packet. The centroid of the packet is at r and at 
k r in wavenumber space. The dispersion relation (4) 
gives uj r = oj r (r,k r ) = const, which has the role of 
a Hamiltonian for the packet. The influence of the 
imaginary part of u> is discussed later. The Hamilto- 
nian equations are 

dr duj r dk r duj r 

~di = V9 = ~dk~ r ' ~df = ~ Ih ' (8) 
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where v g is the group velocity (see for example Bekefi 
1966). In general, the r dependence of ui r results from 
the r— dependence of all of the quantities, ft, k, k critl 
etc., which enter in equation (4). To make the discus- 
sion tractable we assume that Q/k, Q», c s /a r , and £ 9 
are independent of r. Then we have uj r = f2(r)/(fc r ), 
where k r = k r /k cr it{r) and where / is a dimensionlcss 
function. Notice that Figure 4 gives / = w r /0 as a 
function of k r . From equations (8), we find 



dr il df dk r 

dt k crit dk r ' dt 



kcrit ^ dv 



For a flat rotation curve galaxy, Q = £l (r /r) with 
£lo,r = constants, we find dk r /dt = iv r /{k cr itr), and 
r = r (n o /uj r )f(k r ). We may take r to be the ini- 
tial radius of the wave packet. Therefore, for the case 
of Figure 4a, we see that the packet moves radially 
outward in the unstable range of k r and that k r in- 
creases monotonically in that u> r > 0. It is a trailing 
spiral wave. (For k r < 0, f(k r ), which is an even 
function of k r , decreases with k r in the corresponding 
unstable region and the packet moves inward while k r 
increases. The wave in this case is a leading spiral.) 

The maximum amplification of a wave packet re- 
sults when k r increases through all of the unstable 
range of k r in which f(k r ) > (or < 0). This range 
is denoted Ak r . This gives an amplification factor 

A = cxp ( / dt u>i) = cxp ( ^Ll f dk r fccrit ^4r^ 

J J A k r H 

(10) 

In general, A will be less than A max = exp[27r(T 9a ;/ 
T r ot) max], where Tg a i is the age of the galaxy, 
T, r ot is the rotation period at the considered location 
(r), and (u>i/Q) max is the maximum growth rate (as 
a function of k r ). 

For the case m = 2 of Figure 4a, the argument 
of the exponential is about 0.01(k cr itr o )- For values 
of k C rit of the order of that estimated for our galaxy 
w 2n/10kpc at r — 8.5kpc (BT), we conclude that 
the wave growth is insignificant. Basically, the per- 
turbation convects out of the unstable k r range before 
there is appreciable growth. The same conclusion ap- 
plies to the unstable m = 3 mode crossings for ( 9 <1 
and c s /a r = 1. 

3.3. Co- Rotating Stars/High-Mass Counter- 
Rotating Gas 

If the mass fraction of gas £ g is not small com- 
pared with unity, the axisymmetric stability of the 



disk is affected by both the sound speed in the gas 
(c s ) and the radial velocity spread of the stars (ay or 
Q*). For m — 0, equation (4) is a function of uj 2 , and 
the stability/instability boundary uj 2 = gives 



= 1 



l-exp(-X*))/ (X*) 



£>g I 



i + (k r c s / K y ' 

where X* = (k r a r /K) 2 (Toomre 1964). For a given a r 
and £ g , the absence of solutions of equation (9) for k r 
for a critical, sufficiently large c s implies axisymmetric 
stability. This is shown in Figure 5. In contrast with 
the propagating or convective instability discussed 
above, the m = instability gives Wj > with u r = 
and therefore is a non-propagating or absolute insta- 
bility (see for example Lifshitz and Pitaevskii 1981). 
We assume axisymmetric stability. (If this is not the 
case, instability appears in the m > 1 waves giving 
growth independent of the counter-streaming.) 

Figure 6a shows the one-armed (m = 1) star and 
gas mode lines for a high mass fractions of gas and 
c s /a r — 0.316. For the near-crossing shown in Fig- 
ure 6a, the mode lines 'attract' in the vicinity of the 
circle to give instability as shown in Figures 6b and 
6c. (Conversely, the mode lines 'repel' in the case 
where the crossing would be stable.) Figure 7 shows 
the two-armed (m — 2) mode lines and instability 
also for a high mass fraction of gas. Figure 8 shows 
the dependence of the maximum growth rate of the 
one-armed mode on £ g , Q„, and c s ja r - 

Consider the consequence of the wave growth for 
the cases shown in Figures 6 (m = 1) and 7 (m = 2).. 
We comment first on the case of Figure 7 which is sim- 
ilar to the case discussed in the previous sub-section. 
Because u r < 0, f(k r ) < and dk r /dt < over 
the entire unstable range of k r . Over most of this 
range dr/dt > for k r > 0, whereas dr/dt < 
for k r < 0. Thus an unstable trailing spiral wave 
moves radially outward, whereas an unstable leading 
spiral wave moves inward. The maximum amplifica- 
tion from equation (10) is about exp[2.4(fc cr i t r Q )] for 
the trailing waves and exp[1.7(fc cr itr )] for the leading 
waves. Thus, for k cr u = 27r/10fcpc and r Q — lOkpc, 
the amplification factor for the trailing waves is about 
4.2 x 10 6 . 



The case of Figure 6b for m = 1 and 



0.25 



is different because u{k r ) changes sign at k within 
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the unstable range of k r . Recall that for a flat rota- 
tion curve, r — r (Ct/u> r )f(k r ), where r a is the ini- 
tial radius of the wave packet. Thus, for k r < k 
in the unstable range and k r > 0, we must have 
uj r = const. > 0, so that f(k r ) > and dk r /dt > 0. 
However, as k r increases, f(k r ) eventually decreases 
and approaches zero so that r — > 0; that is, the wave 
packet moves towards the center of the galaxy. On 
the other hand, for k r > ko in the unstable range, 
we must have u r — const. < so that f(k r ) < 0, 
and dk r /dt < 0. In this limit, k r decreases and f(k r ) 
eventually increases and approaches zero at k so that 
again r — > 0. For both k r > ko and < k$, the wave 
amplification is even larger than than that for the 
above-mentioned m = 2 instability for the same £ s . 
Thus, it is likely that A = A max . The behavior of 
the trailing waves, k r < 0, is different. For both 
k r < ~k and > — fc , k r moves away from — ko and 
the wave packet moves from its initial radius to larger 
distances. The wave amplification from equation (10) 
in this case is even larger than that for the trailing 
waves because |f2 /ay| ^ 1 in equation (10). When 
the wave-packet reaches and passes the right or left- 
hand limit of the unstable range of k r for k r < 0, the 
packet continues to propagate adiabatically (without 
growth) away from ko with part of the wave energy in 
the upper frequency (w r ) branch and part in the lower 
branch. For example, k r larger than the right-hand 
limit the wave on the lower frequency branch moves 
towards the center of the galaxy whereas the wave on 
the upper branch moves to larger radii. 

For the case of Figure 6c for to — 1 and £ g — 
0.15, uj r is positive throughout the unstable range 
of k r . The maximum amplification factor for the 
trailing waves (k r > 0) is about exp[3.5(fc cr ,t?" )], 
whereas for the leading waves (k r < 0) it is about 
cxp[8.67(fc crit r )]. Both factors are very large for 
kcritfo ~ 2-7T, but that for the leading waves is strongly 
dominant. 

Consider the m = 1 wave growth for large mass- 
fractions of gas, say, £ 9 = 0.3 — 0.7. Stability of the 
to = Toomrc mode requires larger values of 
and/or c s /a r than considered above (see Fig. 5), and 
here we consider Q* = 1.8 and c s /a r = 0.5. For £ g = 
0.3, the unstable waves occur for k r = 0.72 — 1.46, the 
maximum growth rate is Wi/Cl — 0.17, u> r /Q varies 
from 0.092 to 0.096. For £ g = 0.5, the unstable 
range is k r = 0.63 — 2.0, the maximum growth rate is 
Ui/Sl = 0.265, and u r /fl varies from zero to —0.042. 
For £ g = 0.7, the unstable range is k r = 0.64 — 2.56, 



the maximum growth rate is uJi/Q, = 0.29, and u) r /Q 
varies from —0.13 to —0.2. Although the growth rate 
is somewhat larger for £ g = 0.7, the wave amplifica- 
tion is the maximum possible (A max ) for the £ s = 0.5 
case where the wave frequency ui r (k r ) goes to zero in 
the unstable range of k r (see equation 10). 

In Section 5 we discuss non-linear processes which 
may lead to saturation of the wave growth. 

3.4. Co- Rotating Stars/Low-Mass Counter- 
Rotating Stars 

The dispersion relation for a two component galaxy 
consisting of co-rotating (+0) stars and a mass frac- 
tion < i of counter-rotating stars (—SI) is 

= e(u,k r ) = l + (l — £*)T > *(s,k r ) + £*'P*(8 + w,k r ) , 

(12) 

where w = 2mSl/n and "P* is defined below equation 
(2). The co- rotating star modes are given by = 

1 + (1 — ^)V*(s,k r ), while the counter-rotating star 
modes are given by = 1 + £*P*(s + w, k r ). 

Figure 9a shows the behavior of the two-armed 
(to = 2) modes for ^, « 1 and Q» = 1.4 for both 
components. The circled mode crossings in this fig- 
ure satisfy equation (7) suggesting instability. The 
instability is confirmed by Figure 9b which shows the 
complex uj as a function of k r obtained from equation 
(12) for = 0.05 and Q* = 1.4 for both components. 
From equation (10), the largest amplification factor 
is about exp[0.18(fc cr itr o )] for the wave with u> r < 0. 
For k cr i t r <~ 2ir, this amplification is insignificant. 
For the same conditions, the three-armed spiral waves 
are stable. The one-armed waves are unstable for 
k r = 0.935 to 1.5, with a maximum growth rate of 
oj t /Vt = 0.079 at k r = 1.25 where uj r /Q = 0.27, and 
with maximum amplification factor exp[0.13(fc cr ,tr o )], 
which is insignificant for k cr i t r ~ 2tt. 

3.5. Comparable Mass Co-/Counter- Rotating 
Stars 

Figure 10a shows the behavior of the two-armed 
(to = 2) modes for = i and Q» = 1.4 for both 
components. Again the circled mode crossings in this 
figure satisfy equation (7) suggesting instability. The 
instability is confirmed by Figure 10b which shows 
uj r and u>i as a function of k r obtained from equation 
(12). From equation (10), the maximum amplification 
factor is about exp[0.58(fc c ritr o )]. Note that for = 
5, the waves are completely symmetric under 4> 
—<p: the co- rotating wave (w+ = uj r + i^i , w r > 
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0) is matched by an equivalent counter-rotating wave 
(u>- = —Lu r + iu>i) with the same growth rate. Figure 
11 shows the dependence of the maximum growth rate 
on Q*. 

The three-armed waves are stable for = i and 
Q* = 1.4 for both components. However, the one- 
armed waves are strongly unstable for the same con- 
ditions. The mode lines shown in Figure 12a do 
not cross but 'attract' (as noted earlier) near the 
circled point to give instability as shown in Figure 
12b. The exact symmetry <f) — > — for = | 
makes it evident that the co- and counter-rotating 
mode lines in Figure 12a must merge at s — —Vl/n 
to give uv = (independent of Q/k) in the range 
of k r where u>i > 0. Thus the to = 1 instability 
in this case is an absolute instability. In this re- 
spect it is similar to the to = instability. Trail- 
ing and leading wave perturbations can be super- 
posed to give a standing wave in place of equation 
(1), (5$ = Ccos(/ r dr'k r {r')) sin(</>) exp^t). Figure 
13 shows the dependence of the maximum growth rate 
on Q* for = \. 

For = ^, but unequal velocity spreads of the 
components, the symmetry (f> — > —<j> is spoiled. For 
example, the counter-rotating stars may be younger 
with a smaller velocity spread. For u~ = 0.316cr+ 
(Q* = 0.316Q+), we find that stability of the to = 
Toomre mode requires > 1.86, while the to = 1 
mode is unstable for < 2.85. At iQ+ = 1.86, the 
to = 1 mode has uii/il « 0.34 and u) r /Q ~= —0.2. In 
this case the m = 1 instability is convective. 

For < |, the one-armed waves still have large 
growth rates if say > 0.05. This is shown in Figure 
14a. However, uu r takes on positive values due to the 
stronger 'pull' of the co-rotating stars. With ui r > 0, 
k r increases through the unstable range (see equation 
9) for both trailing (k r > 0) and leading (k r < 0) 
waves. Thus the instability is convective. There is 
a weak dependence of u> r /fl on k r which is such that 
the centroid of a wave packet moves outward and then 
inward over a small range of r as k r increases through 
the unstable range. The maximum amplification fac- 
tor from equation (10) is shown in Figure 14b. 

For > i, the solutions are given by those 
for < i by the replacements — > 1 — and 
(uv, u>i) — > (— w r , u>i), provided that the velocity 
dispersions (the Q'^s) of the two components are the 
same. 



3.6. Co- Rotating Stars/Counter- Rotating Gas 
and Stars 

The dispersion relation for a three component 
galaxy consisting of co-rotating stars and a mass frac- 
tion £ g of counter-rotating gas and a fraction of 
counter-rotating stars is 

= e(uj,k r ) = l+{l-£, g -£,*)'P*(s,k r )+£ i g'Pg(s+w,k r ) 

+£*P*(8 + W,kr) , (12) 

where w = 2mfl/K and V* is defined below equation 
(2). The co- rotating star modes are given by w 

1 + (1 — £ 9 — k r ), while the counter-rotating 
gas modes are given by = 1 + £ g Vg(s + w, k r ) and 
the counter-rotating star modes are given by = 1 + 
£*V*(s + w, k r ). 

Figure 15 shows the real and imaginary parts of w 
as a function of k r for one-armed waves for a sample 
case where £ g = = 0.1, with = 1.4 for both star 
components, and c s /a r = 0.316 for the gas. The in- 
stability arises mainly from the interaction of the co- 
rotating stars and the counter-rotating gas due to the 
lower velocity dispersion of the gas. The maximum 
amplification factor for the trailing waves (fc r > 0) 
is about exp[1.23(fc C ritr )], whereas for the leading 
waves (fc r < 0) it is about exp[1.7(fc cr iir )]. For 
kcritr — 27r, these factors are 2.3 x 10 3 and 4.3 x 10 4 , 
respectively. 

Figure 16 shows the dependence of the one-armed 
(to = 1) wave instability on £ g for three values of £*. 
In this figure, Q* = 1.6 for both stellar components, 
and c s /(T r = 0.5. From the discussion of Sec. 3.3, it is 
clear that the largest amplification occurs when uj r = 
within the unstable k r range. Therefore, it follows 
from this Figure 16b that the largest amplification 
occurs for £ ff + k, i. Under this condition the 
amplification is largest for the leading waves (k r < 0). 
The growth rates and amplification factors are smaller 
for the to = 2 waves. 

4. POSITIVE/NEGATIVE MODE ENERGY 
OF UNSTABLE WAVES 

The energy- density of an electrostatic wave in a 
homogeneous plasma is given by the well-known ex- 
pression £ w = uj r (de/duJr)\SE\ 2 /(8ir), where e is the 
dielectric function, (5E = -V<5$ e , and V 2 <5$ e = 
-4ir5p e (Coppi, Rosenbluth, & Sudan 1969). The 
corresponding expression for a self-gravitating medium 
differs by a minus sign and is £ w = — ui r (de/duj r )\Sg\ 2 / 
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(SirG), where Sg = -V£$ and V 2 (5$ = +4irGSp. 
The generalization to tightly wrapped waves in a flat 
galaxy gives the energy (per unit area of the disk) 

£ w = -u r {de/du r )TtG\8Y,\ 2 /\k r \ , (13) 

where we have used <5<f> = —2-wG5Y,/\k r \, and 5T, is the 
total surface mass-density perturbation. An equation 
equivalent to (13) is derived by Shu (1992). 

For a galaxy consisting of co- and counter-rotating 
components (now denoted for generality by + and — ), 
we have e = 1 + V + + V~ . Thus, the wave energy- 
density is the sum of the energy densities associated 
with the modes in the co- and counter-rotating com- 
ponents, 

f — f + + f~ 

£± = -u r {dV ± /du r )TrG\5Y,\ 2 /\k r \ . (14) 

If the signs of £+ and £~ arc different, we have 
the necessary condition for the well-known instabil- 
ity of interacting positive and negative energy modes 
(Coppi et al. 1969). The negative energy mode can 
grow by feeding energy into the positive energy mode. 
A plasma example is the two-stream instability which 
can occur when a beam of charged particles passes 
with velocity Vb through a background plasma. In 
this case the energies associated with the modes in 
the beam and the background plasma have opposite 
signs. The mechanism of this instability can be un- 
derstood by considering a localized excess or clump of 
say positive charge in the beam which will induce a 
negative clump in the background. The relative mo- 
tion of the clumps leads to an electrostatic attraction 
between them, a slowing down of the beam clump, a 
speeding up of the background clump, and a growth 
of the electric field energy. The clump size must be 
sufficiently large to give instability, larger than about 
■KVb/uip for the case of equal beam and plasma densi- 
ties, where u> p is the plasma frequency In a homoge- 
neous self-gravitating system involving a beam and a 
background, a density excess in the beam will induce 
an excess in the background, and it would appear that 
a similar two-stream instability would occur. How- 
ever, for most distribution functions, the large clump 
size needed for the two-stream instability (> nvb/wj, 
where ujj = (AirGp) 2 is the Jeans angular frequency) 
results instead in the Jeans instability being dominant 
(Lynden-Bell 1967, Araki 1987). The Jeans instabil- 
ity occurs for clump sizes larger than about ttAv/ujj, 
where Av is the velocity spread. In contrast, in a disk 
galaxy the velocity spreads of the ±fi components and 



the disk rotation [(k + ) 2 = (k ) 2 > 0] provide stabil- 
ity to the Jeans instability (m = 0) (Toomre 1964). 
From equations (2) and (3), the condition for 
< is simply s r (s r + w) < or — O < uo r /m < 
ft which is equation (7). Recall that u r /m is the pat- 
tern or angular phase velocity of the wave. For the 
two-stream instability in a beam-plasma system, the 
phase velocity of the unstable waves is also interme- 
diate between velocity of the beam and that of the 
background. 

The axial angular momentum of the wave (per unit 
area of the disk) can be written as J w = + J~ , 
where J± = -m{&P ± /dw r )7rG|(5£| 2 /\k r \ (Coppi et 
al. 1969), so that £ w = {uj r jm)J w . Thus, the wave 
angular momenta associated with the modes in the co- 
and counter- rotating components have opposite signs 
under the same condition that the mode energies have 
opposite signs. Under this condition, the angular mo- 
mentum of the mode in the co-rotating component is 
negative, whereas that in the counter-rotating compo- 
nent is positive. Thus, the total angular momentum 
of the mode plus matter of the +fi component is re- 
duced, whereas that in the — Q component is increased 
(but decreased in magnitude). The situation is analo- 
gous for the beam-plasma instability where the linear 
momentum of the beam mode is negative and that of 
the background mode is positive. 

5. NON-LINEAR EFFECTS AND GAS AC- 
CRETION 

The non-axisymmetric instabilities discussed here 
are of interest because they could be important for 
the inward radial transport or accretion of counter- 
rotating (— ft) matter deposited or accreted at large 
r onto an existing flat galaxy consisting mainly of co- 
rotating (+ST) stars. We consider first the case where 
the counter-rotating matter is gas. 

The linear wave growth is determined mainly by 
three parameters, the Toomre parameter for the 
stars, the ratio of the sound speed in the gas to the ra- 
dial velocity spread of the stars c s /a r , and the mass 
fraction of the counter- rotating gas £ g . In general 
these quantities will vary with r; the theory is still 
valid as long as the variation is small on a scale l/k r . 
For a given £ g , the wave growth rate decreases mono- 
tonically with increasing and c s /a r , as shown for 
example by Figure 8 for one-armed waves. On the 
other hand, with Q* and c s ja r fixed, the wave growth 
is zero for £ s below a threshold value and increases 



8 



strongly as £ s increases above this value as shown by 
Figure 8a. 

Suppose that Q*, c s /oy, and £ g are such as to give 
a large growth rate, say, Wj/Q > 0.05, and an appre- 
ciable wave amplification factor, A^> 1. After several 
rotation periods of the matter at a radius r, the wave 
amplitude can grow to a level where non-linear effects 
become important (|<5$| - (0.01 - 0.05)(fh-) 2 ). One 
non-linear effect is the scattering of stars into less cir- 
cular orbits which acts to increase the radial velocity 
spread and therefore Q*. In turn, the increase of Q* 
acts to reduce the growth rates. However, a finite level 
of non-axisymmctric waves may remain excited due 
to reduced but non-zero wave growth. These waves 
can give an effective viscosity v e which causes inward 
transport or accretion of the gas and outward trans- 
port of the gas's angular momentum. The gas can 
remain in approximately circular motion by radiat- 
ing the energy excess which results from its inward 
motion. This gives an accretion luminosity (per unit 
area of the disk) T, g Q 2 r\v gr \, assuming a flat rotation 
curve, where v gr < is the average accretion speed of 
the gas. A rough estimation using quasi-linear theory 
gives v e ~ Q,r 2 \k r r\ 3 \S^>/(flr) 2 \ 2 and v gr ~ —v e /r, 
where <5$ is the residual wave level. Relevant values 
for accretion from say r = 20 kpc in 3 x 10 9 yrs. are 
\v gr \ > 5 km/s. 

A second non-linear effect is the steepening of the 
spiral wave in the gas which can lead to the forma- 
tion of a oblique shock as indicated in Figure 17 for 
the case of a one-armed leading wave. (The amplifi- 
cation factor for the leading wave is generally larger 
than that for the trailing wave.) The gas velocity 
parallel to the shock is unchanged across the shock. 
Up stream of the shock, the normal component of 
the velocity v gn = 9£lr is larger than the ambient 
sounds speed c° in order to have a shock. Here, 
= |/c r r| _1 <C 1 is the pitch angle of the spiral arm. 
That is, the normal Mach number M n = v gn /c° > 1. 
We assume that M n is not much larger than unity 
so that the shock is not very strong. Downstream 
from the shock, the normal velocity of the gas is less 
than the post-shock sound speed. As a result, the 
gas acquires an inward radial velocity v gr = — /3f2r, 
where (3 = — |^-(1 — 1/M 2 )9 and 7 is the adiabatic 
index. This shock deflection could be important for 
the inward radial transport of counter-rotating gas. 
The angular momentum lost by the gas may in be 
transported outward by the spiral wave. Also in this 
case, the gas can remain in approximately circular 



motion by radiating the energy excess arising from ac- 
cretion. Across the disk, the azimuthally averaged ac- 
cretion speed v gr has a strong non-linear dependence 
on £ g because at small t; g there is no instability and 
no transport, while for £ g m \ there is maximum wave 
amplification (Sec. 3.3). As a result, the gas trans- 
port equation dT, g /dt + (l/r)d[rT, g v gr (T, g ))/dr = 
(neglecting star formation) is strongly non-linear. Its 
solutions may involve inward propagating soliton-like 
perturbations in T, g . This may be the cause of ring- 
like features observed in some counter-rotating galax- 
ies (for example, NGC 7217, Buta ct al. 1995, and 
NGC 4138, Jore et al. 1996). 

The growth of spiral arms with significant self- 
gravity may be important in inducing collisions and 
agglomeration of gas clouds in the arms (see for exam- 
ple Roberts, Lowe, & Adler 1990) thereby enhancing 
the rate of formation of counter-rotating stars. Newly 
formed counter-rotating stars should have a smaller 
velocity spread than the older co-rotating stars which 
have undergone stochastic heating (BT, p. 484). 

If the counter-rotating matter consists of stars, the 
m = 1 instability discussed in Sec. 3.2 may be im- 
portant. From Figure 14, the instability is significant 
for a mass fraction of counter-rotating stars larger 

than 0.05 and is strongest for — \ if Q* = 1.4 

^ 1 

for both components. For increasing from ^, the 
growth rate decreases. The large growth rates and 
amplification factors of Figure 14 we expect to ex- 
ist only transiently because the resulting waves would 
rapidly scatter stars into less circular orbits increasing 
Q* and decreasing uji. However, a finite level of waves 
may remain excited due to much smaller growth rates. 
Consider now the possibility of radial inflow of the 
stars. A counter-rotating star at a radius r with an- 
gular momentum l~ < can lose angular momentum 
Sl~ > to the spiral wave in a given period of time 
so that the magnitude of l~ decreases. At the same 
time, a co-rotating star also at r with angular momen- 
tum > can lose <5Z+ < so that the magnitude of 
If also decreases. The resulting smaller values of \lf\ 
correspond to smaller mean radii. In order for a co- or 
counter-rotating star to stay in approximately circu- 
lar motion, its energy must decrease by an amount 
Se^ — Q5lf < 0, assuming a flat rotation curve. 
However, the ratio of the energy to the angular mo- 
mentum of the spiral wave E w j J w — u r / m = uo r <C f2 
from Sec. 4 and Figure 14. That is, the energy loss 
of the stars to the spiral wave is too small to allow 
inward radial motion on approximately circular or- 
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bits. Even a small fractional inward motion of the 
stars will increase sufficiently to stabilize the spi- 
ral wave modes. 

The stability of a galaxy with both counter-rotating 
stars and gas (Sec. 3.6) appears qualitatively similar 
to the results found with only counter-rotating gas. 
The largest wave amplifications are for the m = 1 
leading spiral waves and for + £ g w \ . The above 
discussion of gas transport is also pertinent to this 
case. 

6. CONCLUSIONS 

The two-stream instability - well-known in plasma 
physics - has been discussed by a number of authors 
for counter-streaming self-gravitating systems. How- 
ever, a somewhat muddled picture has emerged for 
the role of this instability in galaxies with counter- 
rotating components including the notion that the 
instability does not occur due to the difference be- 
tween the universal attraction of masses and the re- 
pulsion/attraction of like/unlike charges. Indeed, sta- 
bility analysis of a homogeneous self-gravitating sys- 
tem of counter-streaming particles (following Jeans' 
neglect of the background potential) , shows that the 
two-stream instability is dominated by the Jeans in- 
stability if the distribution functions are Maxwellian 
(Lynden-Bell 1967, Araki 1987). However, the homo- 
geneous system results are not relevant to counter- 
stream matter in disk galaxies where the Jeans in- 
stability (m = 0) is stabilized by both the veloc- 
ity spreads of the two components and the disk 
rotation (k 2 > for both components) (Toomre 
1964). A study by Araki (1987) of the stability 
of counter-streaming, rigidly-rotating, finite-radius 
Kalnajs disks docs in fact show anm=l two-stream 
instability, but the physical nature of the instability is 
not elucidated. Simulation studies of different galaxy 
models by Zang and Hohl (1978), Sellwood and Mer- 
ritt (1994), and Howard et al. (1996) all show clear 
evidence of a strong m = 1 two-stream instability. 

In this work, we first briefly summarize known 
results for small-amplitude, tightly-wrapped spiral 
waves in a single component stellar disk and show (in 
Sec. 2) that the higher order modes of the dispersion 
relation have been overlooked in standard treatments 
(BT, Palmer 1994). The dispersion relation gives the 
dependence of the wave frequency (co) on the radial 
wavenumber (k r ), and the omitted curves are labeled 
I = 2,3,.. in Figure 2. These higher order modes 



are the analogues of the Bernstein modes in a col- 
lisionless plasma which propagate across a uniform 
magnetic field. Their existence leads to a richer set 
of Linblad resonances which occur at pattern speeds 
Cl p = uj/m = 0±//t/m with m = 1,2, .. and I — 1, 2, ... 

We assume tightly wrapped spiral waves, \k r \r » 
1, so that the WKB approximation holds, but it is 
not known how strong this inequality must be. On the 
other hand, for effective swing amplification involving 
radial reflection of waves, which is thought to drive 
spiral waves in galaxies with only co-rotating matter, 
one needs \k r \r < 3 (BT, p. 376). Here, the two- 
stream wave amplification is sufficiently large so as 
to be important without wave reflection. Although 
selection effects may be involved, observed counter- 
rotating disk galaxies are early-type spirals or SO's 
which do not show prominent open spiral structure. 

We go on to discuss spiral waves in flat galaxies 
with counter-streaming where there is a co-rotating 
stellar component with angular velocity +CI and a 
counter- rotating component (s) of gas and/or stars 
with angular velocity —Cl. We then plot on the 
same graph the oj = uj(k r ) lines - referred to as 
mode lines - for the co-rotating stellar component (ne- 
glecting the counter-component) and the w = cj(k r ) 
lines for the counter-component (neglecting the co- 
component). At crossing points of these mode lines 
in the (k r , u) plane, a strong resonant interaction can 
occur between co- and counter-rotating components. 
Direct solution of the dispersion relations (Sec. 3) 
shows that the crossings are unstable if the pattern 
angular speed of the spiral wave u>/m is such that 
— Cl < u/m < CI. The 'mode-crossings' and near 
'mode-crossings' (Sec. 3.3) satisfying this condition 
give the two-stream instability in a counter-rotating 
galaxy. They arise from the fact that the mode en- 
ergies for the co- and counter-rotating components 
have opposite signs for —Cl < uj/m < Cl (Sec. 4). 
This is analogous to the two-stream instability in a 
beam-plasma system where the mode energies for the 
beam and the background have opposite signs for a 
wave phase velocity intermediate between that of the 
beam and that of the background. 

Instability of a spiral wave may or may not be sig- 
nificant because the waves are in most cases convec- 
tive in the sense that both the position (r) and the 
wavenumber (k r ) of an unstable wave packet evolve 
with time (Sec. 3.2). More important than the lo- 
cal growth rate (wj) is the amplification factor of the 
advected wave A (equation 10) which is the maxi- 
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mum factor by which an unstable packet can grow 
before being convected out of the unstable range of 
k r . In most cases, we find that the amplification fac- 
tors are largest for one-armed m — 1 leading spiral 
waves (with respect to the co- rotating stars). How- 
ever, for larger stellar velocity spreads (larger Q*), 
the two-armed m = 2 wave may be unstable, with 
the largest amplification for the trailing wave, while 
the m = 1 wave is stable. The growth rates and 
amplification factors increase as the mass-fraction of 
counter-rotating gas £ g or of stars increases if these 
fractions are less than ss i. On the other hand, both 
u>i and A decrease as the velocity spreads of the stars 
and/or gas increase. 

For the case of only counter-rotating gas of mass 
fraction £ g <« |, the largest amplification factors 
are for the m = 1 leading spiral waves (Sec. 3.3). 
For the case of only counter-rotating stars of mass- 
fraction <« \, the largest u>i and A values occur 
also for the m = 1 leading spiral waves (Sec. 3.5). 
If the velocity spreads in the two stellar components 
are equal, then for = |, the m = 1 instability is an 
absolute or non-propagating instability, the real part 
of the wave frequency (u r ) is zero, and the growth rate 
has its maximum value. For > ^, the instability is 
again convective, the amplification factors are largest 
for the trailing spiral waves, and u>i and A decrease 
with increasing £*. These results are in qualitative 
accord with the theoretical results of Araki (1987) for 
counter-rotating stellar disks and with the simulation 
results of Sellwood and Merritt (1994), where = \ 
in both studies. 

For a galaxy with both counter-rotating stars and 
gas, the largest amplification is for the m = 1 leading 
spiral waves, and it occurs when the counter-rotating 
mass fraction is £* + £ g w \ (Sec. 3.6). 

Possible non-linear effects which act to limit the 
wave growth and amplification are discussed in Sec. 
5. The effects include scattering of star orbits by 
the wave which increases Q*, and heating of the gas. 
A residual level of spiral waves may remain excited 
which give an effective viscosity for the gas causing 
its accretion. Also, a leading m = 1 spiral shock 
wave may form in the gas causing its accretion. How- 
ever, from the considerations of Sec. 5, it appears un- 
likely that the spiral waves cause accretion of counter- 
rotating stars. 

The gas viscosity (and thus accretion rate) due to 
the two-stream waves is plausibly an increasing func- 
tion of the spiral wave amplification factor. With 



no amplification there are no waves and no viscos- 
ity. Thus, the fact that the largest amplification 
occurs for a mass-fraction of counter-rotating mat- 
ter + £ g i may be pertinent to the counter- 
rotating galaxy NGC-4550 (Rubin et al. 1992) which 
is remarkably symmetric between co- and counter- 
rotating stars, w \ (Rix et al. 1992). A schematic 
scenario is that initially counter-rotating gas is sup- 
plied at large r and has a fastest accretion rate for 
£ g = ^. At later times, after counter-rotating stars 
have formed with mass fraction £*, the fastest accre- 
tion rate is for £ g = | — £». At even later times, the 
gas accretion ceases when w i. 

Accretion of gas due to two-stream waves is strongly 
nonlinear when £ s is large and this may lead to a 'pile 
up' of gas into rings at one or more radial distances. 
Prominent rings are seen in NGC 3593, NGC 4138, 
and NGC 7217. Buta et al.(1995) have suggested that 
the locations of the three rings in NGC 7217 corre- 
spond to different Linblad resonances. At sufficiently 
high counter-rotating gas densities counter-rotating 
star formation may be triggered. 
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Fig. 1. — Sketch of the counter rotating disk geome- 
try. 

Fig. 2. — The figure shows the frequency [s = 
(lo — mfl) / k] - wavenumber (fc r ) dependence of the 
tightly wrapped spiral wave modes in a disk of stars. 
Here, ft is the angular velocity of the stars, k is their 
epicyclic frequency, k cr u is the critical wavenumber, 
and I labels different branches of the dispersion rela- 
tion. The higher branches (I > 2) are analogues of the 
Bernstein modes in a magnetized plasma as discussed 
in the text. 

Fig. 3. — The figure shows the different star and gas 
two-armed (m = 2) modes in a galaxy of co-rotating 
stars with = 1.4 and a small mass fraction of 
counter-rotating gas (£ s <C 1). The curves extend to 
negative k r as even functions. At the circled points 
where the star and gas modes cross there is a strong 
resonant interaction which may lead to instability. 
Here, only the point B is unstable. The detailed be- 
havior of the mode crossings is shown in Figure 4. 

Fig. 4. — The figure shows the behavior of the 
star/gas mode crossings for two armed waves (m = 2) 
and a small mass fraction of gas, £ g = 0.01. Also, 
= 1.4, c,,/oy = 1, and we have assumed a flat ro- 
tation curve so that k — >/2tt. a is for an unstable 
case (point B of Figure 3), and b for a stable case 
(point C of Figure 3). 

Fig. 5. — The figure shows the dependence of the ax- 
isymmetric (m = 0) Toomre (1964) stability thresh- 
old on the sound speed in the gas (c s ) obtained from 
equation (9). Here, ay is the radial velocity disper- 
sion of the stars, £ g is the mass fraction of the gas, 
and Q* is the Toomre stability parameter for the 
stars defined below equation (2). We have assumed 
k = y/2Q. Note that the dashed curve £ g = 1 is 
given by c s /a r — 0.9355/Q*, which is the same as 
c s = n/(2k cr it) from equations (2) and (3). 

Fig. 6. — The figure shows the nature of the one- 
armed (m = 1) modes in a galaxy of co-rotating stars 
and an appreciable mass fractions of counter- rotating 
gas. For this figure, Q» = 1.4, c s /a r = 0.316, and 
k = \/2Q. a shows the star and gas mode lines. The 
circle indicates the near crossing which is unstable, 
b shows the real and imaginary parts of u obtained 
from equation (4) for £ g — 0.25, and c shows the same 
quantities for £ g = 0.15. 



Fig. 7. — The figure shows the nature of the two- 
armed (in = 2) modes in a galaxy of co-rotating 
stars and an appreciable mass fraction (£ s = 0.25) 
of counter- rotating gas. For this figure, Q* = 1.4, 
C s/°V = 0.316, and n — \/2tt. a shows the star and 
gas mode lines. The circles indicate unstable mode 
crossings, b shows the real and imaginary parts of 
u) obtained from equation (4). 

Fig. 8. — The figure shows the dependences of the 
maximum growth rate of the one-armed waves on £ g , 
Q*, and c s /a r . The left-hand limit in b at c s j a r = 0.2 
is close to the stability threshold of the axisymmetric 
instability (see Fig. 5). 

Fig. 9. — The figure shows the nature of the two- 
armed (m = 2) modes in a galaxy of co-rotating 
stars (+f2) and a small fraction (£* <§; 1) of counter- 
rotating stars (—0). a shows the mode lines with 
the unstable crossings circled, b shows the real and 
imaginary parts of ui obtained from equation (12). 
For both star components, Q* = 1.4, and we have 
taken k = \/2Q. The vertical arrows on the u>i curves 
indicate the associated uy curves. 

Fig. 10. — The figure shows the nature of the two- 
armed (m = 2) modes in a galaxy of co-rotating stars 
and an equal mass (£* = \) of counter-rotating 
stars (— Q). a shows the mode lines with the unsta- 
ble crossing circled, b shows the real and imaginary 
parts of the wave frequency ui as a function of k r ob- 
tained from equation (12). For both star components, 
= 1.4. Also, we have taken k = \/2Q. The vertical 
arrows on the un curves point towards the associated 
uo r curves. 

Fig. 11. — The figure shows the Q* dependence of 
the maximum growth rate of the two-armed waves for 
= i. The value of k r for the maximum growth is 
also indicated. We have taken n — 

Fig. 12. — The figure shows the nature of the one- 
armed (m = 2) modes in a galaxy of co-rotating stars 
and an equal mass (£* = |) of counter-rotating 
stars (—0). a shows the mode lines, and the cir- 
cle indicates the unstable near crossing of mode lines, 
b shows the real and imaginary parts of the wave fre- 
quency uj as a function of k r obtained from equation 
(12). For both star components, Q» = 1.4. 

Fig. 13. — The figure shows the dependence of 
the maximum growth rate of the one-armed waves 
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for a galaxy of equal mass co- and counter-rotating 
components, = |. 

Fig. 14. — The figure shows the maximum growth 
rate and amplification factor of the one-armed waves 
for a galaxy of co-rotating stars and an appreciable 
mass-fraction of counter-rotating stars. 

Fig. 15. — The figure shows the nature of the one 
armed unstable mode of galaxy of co-rotating stars 
and a fraction £ g = 0.1 of counter- rotating gas and a 
fraction = 0.1 of counter-rotating stars. For this 
figure, <5„ = 1.4 for both star components, c s ja r = 
0.316, and n = \/20. 

Fig. 16. — The figure shows the dependence of the 
one-armed (m = 1) spiral wave instability on £ g for 
three values of £». a shows the dependence of 
(uj r ) m on £ g , where (to r )m is the value of u> r at the 
point where |w r (fc r )| is a minimum for k r within the 
unstable range. b shows the dependence of the 
maximum growth on The figure assumes Q* = 1.6 
for both stellar components, c s /a r — 0.5, and k = 
\/2Q. The largest wave amplification occurs when 
(ujr) m = 0, and panel a shows that this happens for 
+ £g ~ 5 • 

Fig. 17. — Geometry of a one-armed leading spiral 
shock wave in counter-rotating gas which leads to in- 
ward radial motion of the gas. 
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